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We study the dynamics of a single Frenkel exciton in a disordered molecular chain. The coherent- 
potential approximation (CPA) is applied to the situation when the single-molecule excitation en- 
ergies as well as the transition dipole moments, both their absolute values and orientations, are 
random. Such model is believed to be relevant for the description of the linear optical properties 
of one-dimensional J aggregates. We calculate the exciton density of states, the linear absorption 
spectra and the exciton coherence length which reveals itself in the linear optics. A detailed analysis 
of the low-disorder limit of the theory is presented. In particular, we derive asymptotic formulas 
relating the absorption linewidth and the exciton coherence length to the strength of disorder. Such 
expressions account simultaneously for all the above types of disorder and reduce to well-established 
form when no disorder in the transition dipoles is present. The theory is applied to the case of 
purely orientational disorder and is shown to agree well with exact numerical diagonalization. 
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I. INTRODUCTION 

The optically active states in organic molecular ag- 
gregates^ — the spatially-regular linear arrangements 
of dye molecules — are one-dimensional Frenkel exci- 
tons£ The nature of the excitonic states to be extended 
over many molecules reveals itself in the narrowing of 
the absorption resonances and the shortening of the ra- 
diative lifetime upon aggregate formation (effects known 
as exchange narrowing and "superr adiance")^ The re- 
lated coherence length coincides with the total number of 
molecules in the aggregate only if the size of the latter is 
sufficiently small. Otherwise various exciton decoherence 
mechanisms, among which the static disorder is the most 
essential at low temperature, substantially reduce the co- 
herence length that becomes independent from the actual 
size of the aggregate. Considerable work has been done 
to analyze the aspects of disorder in the molecular ex- 
citation energieo 4 i 5 i 6 i 7 i 8 i 9 i 10 i 1:L and positionsi 12 i 13 i 14 i 15 i 16 
Much less attention has been devoted to the disorder in 
orientations^ or, generally, in the transition dipole mo- 
ments of individual molecules. However, from the limited 
available information on the actual structure of molecu- 
lar arrangements one cannot exclude that the transition 
dipoles are seriously influenced by disorder. Moreover, 
the last type of disorder affects in a nontrivial way the 
separate components of the optical susceptibility tensor 
and, hence, would be observable in polarization-resolved 
absorption and luminescence experiments 18 . The present 
work has been motivated by the need for a detailed the- 
oretical analysis of orientational disorder, as well as clar- 
ification of certain aspects concerning the exciton coher- 
ence length in molecular aggregates. 

In this paper the linear optical properties of a disor- 
dered molecular chain are studied using the single-site 



coherent-potential approximation 19 - 20 (CPA). Applied to 
various disorder problems, such as that of the electronic 
structure in a random alloy, this self-consistent approx- 
imation was shown to reproduce well the single-particle 
characteristics like the density of states (DOS). As for 
the systems in which elementary excitations are Frenkel 
excitons, the CPA is in addition capable to provide the 
complete information on the linear optical response, be- 
cause the latter is extracted only from the single-exciton 
Green's function (GF) averaged over disorder realiza- 
tions. The CPA has been successfully used to model the 
optical spectra of periodic molecular arrangements with 
random on-site energies 21 * 2 ^ and composition. 23 How- 
ever, if the transition dipoles are also affected by dis- 
order, a proper modification of the theory is required. 
The key observation which allows us to do this is that 
the dipoles enter the off-diagonal part of the Hamilto- 
nian in a bilinear form. In this case a single-site approx- 
imation can be constructed according to a vector analog 
of the Shiba ansatz^ The resulting scheme is equiva- 
lent to the matrix extension of CPA derived by Black- 
man, Esterling and Berk 25 (BEB) for compositionally- 
disordered alloys with random hopping energies. [The 
analytical proper ties of BEB CPA have been examined 
in Refs. 126.27. 28 29] . Persson and LiebscbJfi, and Rozen- 
baum et al£L have constructed a similar version of CPA 
to study, respectively, the susceptibility of polarizable 
particles and the vibrational modes of coupled oscilla- 
tors randomly and isotropically oriented in two or three 
dimensions. However, their theory is still not applica- 
ble to the molecular aggregates in which one typically 
deals with an essentially anisotropic distribution of the 
transition dipoles. 

The paper is organized as follows. In the next section 
we define the model of a disordered molecular aggregate 
and set up the basic formalism for the forthcoming cal- 
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culations. The exciton Green's tensor is introduced to 
express the related single-particle quantities: the DOS, 
the linear polarizability, and the coherence length. In 
Sec. IIIII we present the general formulation of the BEB- 
CPA scheme and in Sec. IIVI derive its analytic solution 
in the low-disorder limit. In Sec. E] the theory is applied 
to a specific case of orientational disorder and compared 
to the results of the exact numerical diagonalization. 



II. BASIC FORMALISM 
A. Model 

An aggregate can be considered as a chain of Af iden- 
tical two-level molecules. The linear optical properties in 
the resonance region are determined by the single-particle 
states of the Frenkel Hamiltonian^ 



H 
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The operator B^ (B n ) creates (annihilates) an electronic 
excitation of energy e n on the nth molecule (whose ex- 
cited state is supposed to be nondegenerate). The trans- 
fer terms J nm result from the Coulomb interaction be- 
tween an excited molecule and one in the ground state. 
Accounting only for the dipole-dipole contribution one 
has 
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where p" are the vector components of the nth molecule 
transition dipole. The "coupling kernel" i^nm carries in- 
formation on the dielectric function of the surrounding 
material and the location of the molecules, but not on 
their dipole moments. If the dynamics of the exciton sub- 
system is nondissipative (and the molecules have no mag- 
netic structure, so that p" can be chosen real) "d"^ has 
to be symmetric in the tensorial indices {d^m — $nm)- 
We shall assume the molecules to form a one-dimensional 
regular lattice such that depends solely on the in- 
termolecular distance |n — m\ (measured in the units of 
lattice spacing). The theory will be presented in a gen- 
eral way without specifying a concrete form of this de- 
pendence, i.e., for any screening and anisotropy of the 
Coulomb interaction that can occur in a solvent or on a 
substrate interface. 

The disorder enters in our model both through the 
on-site energies e„ and the dipole moments p". We shall 
assume the families of random parameters {£ n ,Pn} corre- 
sponding to distinct sites to be mutually independent and 
to have identical probability distribution. At the same 
time we still allow for correlations between the transition 
energy and the dipole components of a given molecule. 

Finally, we shall neglect everywhere the effects of a 
finite length of the aggregate considering the limit M — > 
oo. This approximation is justified for sufficiently strong 



disorder when the exciton coherence length found for the 
infinite system does not exceed the actual length of the 
chain. 



B. Exciton Green's tensor 

For a given random realization the total information on 
the aggregate's optical dynamics is imbedded in the com- 
plete set of single-exciton eigenstates J2 n ipn(E)B} l |0); 
the vacuum |0) is the direct product of the molecular 
ground states; ip n {E) is the real-space exciton wave func- 
tion with an eigenenergy E. The local density of states 
(LDOS) on the nth site and the total (normalized) exci- 
ton DOS are found as 

p n (w) = V \ME)\ 2 S(cu — E) = -~ImG nn (w), (3) 

{E} 



{E} 
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where the real-space single-exciton GF is given by 
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and the symbol {E} indicates that summation is per- 
formed over all eigenstates. Here and below the GF for a 
real- valued argument lu is found as the limit z — > uj + iO + . 
Assuming that around the single-exciton resonance the 
light wavelength is much larger than the spatial extent 
of a typical excitonic wave function the aggregate inter- 
acts with the optical field as a pointlike dipole. Hence, 
the linear response to the external optical field is com- 
pletely described by the polarizability tensor 

r+<x> 

X aP (z) = i dte izt (0\[P a (t),P^(0)}\0). (6) 
Jo 

in which the operator of total polarization represents the 

sum of local polarizations, P a = J^ni^'"^ + Pn )> 
both positive- and negative-energy components: 



pa(+) _ „<*Rt p«(-) — n a n 
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To set up a convenient formalism for the forthcoming 
discussion let us define the time-ordered two-point corre- 
lator of the local polarizations 
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Tf m (z) = - l dte^(0\P^-\t)P^ + \0)\0), (8) 
Jo 



or, equivalcntly, 



r"m( Z ) — Pn^nm(z)p^ n 



(9) 
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The newly introduced quantity will be referred through 
the paper as the Green's tensor (GT), as opposed to the 
"scalar" GF considered above. Combining the previous 
formulas we arrive at the straightforward relation 



E 



[Ki(z) + r^*(-z*)], (io) 



which indicates that (up to the sign factor) the GT co- 
incides with the positive-energy counterpart of the local 
linear response function. 

As usual one is interested not in the solution of the 
entire dynamical problem for a given disorder realiza- 
tion, but in finding the configurational average of some 
basic observable parameters. In our case the quanti- 
ties of interest are those defined in Eqs. and (|10fl . 
The subsequent implementation of CPA suggests that 
the disorder-averaged DOS p(u>) = (/o(w)) and polariz- 
ability x a ^( z ) = (x a ^i z )) are to be expressed in terms 
of the conditionally-averaged on-site GT 



all sites except nth 
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with fixed nth site variables {e n ,p"}, and the complete 
statistical average 



(12) 



the only quantities accessible within a single-site effec- 
tive theory. To proceed with the DOS it is sufficient 
to notice that, because the nth-site dipoles p" remain 
fixed, the conditional average of the local GT and GF 
are still related to each other as established in Eq. 
^nn( z ) = PZGnn{z)pP- As a result the scalar GF is ob- 
tained via projection 

(13) 

Using this relationship the conditionally-averaged LDOS 
can be expressed as 



\pr, 



(14) 



while the total DOS is found in terms of f"^(z) as a 
trivial single-site average 
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Remarkably, the statistical averaging in the last formula 
can be performed immediately if the probability distri- 
bution of the transition dipoles is of the "purely ori- 
cntational" form. The trace over vector indices of the 
conditionally-averaged GT depends only on the deter- 
ministic, in this case, absolute value of the nth-site dipole 



moment \p n \- The averaged GT and GF are related to 
each other as 

J2Kn^) = \Pn\ 2 G nn (z), (16) 
oc 

and consequently, 

pM = -^^E Imf ««H- 

' 1 a 

The straightforward application of the averaging pro- 
cedure to Eq. IjlOl) gives the disorder-averaged linear po- 
larizability in the form 

x at3 (z) = -M^ KM + (-**)] ■ (is) 

m 

In order to reduce the double real-space summation we 
have used the translational invariancc of T^(z) that al- 
lows to replace in both terms of Eq. (|T%)l J2 n m ^nm( z ) by 
•^Em^nrnf 2 )- The last sum here is already converging 
and does not depend on its upper limit J\f. Equivalently, 



r p {z) = -N r^ (z) + r^(-z*) 



=,a/3 



a/3* 



(19) 



where the momentum-domain disorder-averaged GT is 
introduced according to 



Ki(z) l_ 7r 27T e 



dk 



ik(n-rn) pa/3 



(20) 



k being measured in units of the inverse lattice spac- 
ing. Noticeably, the above polarizability scales linearly 
with the number Af of molecules constituting the aggre- 
gate. Unlike the physics coming from the nontrivial de- 
pendence of Tnm( z ) 011 t ne energy variable z and the 
intersite separation \n — m\ this elementary aspect is not 
related to the degree of exciton coherence. Nevertheless 
it guarantees fulfillment of an important part of the gen- 
eral sum rule according to which the polarizability is an 
extensive quantity proportional to the total number of 
polarizable objects. 



C. Coherence length 

In the literature dealing with the optics of molecular 
aggregates it is usually introduced the notion of the exci- 
ton coherence lengthii*i2ii^ii to characterize the spatial 
extension of the exciton wave functions. The cited works, 
addressing the problem mainly with numerical methods, 
employ several definitions of this parameter to be ex- 
tracted from the set of wave functions obtained with ex- 
plicit diagonalization procedures. Even though the quan- 
titative estimates provided within all approaches can be 
in reasonable agreement with each other, the physical 
arguments used are somewhat different. From the view- 
point of the present work a natural definition of the exci- 
ton coherence length can be given on the basis of Eq. i|18[l • 
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In fact, entering it the disorder- averaged single-exciton 
GT contains all the information on the linear optical re- 
sponse of the excitonic system and depends nontrivially 
on the amplitude and phase coherence between the wave 
functions for different disorder realizations. These co- 
herence properties partly reveal themselves in the DOS 
and absorption spectra which can be essentially different 
from those of the disorder-free system, but it is also the 
dependence of T"^(z) on the intersite separation that 
gets strongly affected by the structural disorder within 
the chain. Excluding some peculiar cases ^ each compo- 
nent of the disorder-averaged GT is characterized by an 
exponential behavior 

f^(z)~exp( l r /3 W|n-m|), \n - m\ -> oo, (21) 

where £, al3 (z) is a complex- valued wave number. Keep- 
ing in mind this asymptotics while performing the real- 
space summation in Eq. (|18|l one concludes that the dom- 
inant contribution to the disorder-averaged polarizability 
X a ^( z ) comes from the pairs of sites with Im£ a ^ (z)\n — 
to | < 1. The parameter 



indicating how many molecules contribute to a given 
component of the linear polarizability, provides a nat- 
ural definition of the exciton coherence length. Clearly, 
N a @(z) is a straightforward generalization of the con- 
ventional quasiparticle phase coherence length, defined 
in terms of the large-intersite-separation asymptotics of 
the GF, to incorporate the vectorial nature of the exci- 
tonic polarization. The quantity 12 21) can be also thought 
of as the nonlocality range of the linear optical response 
function. 



III. BEB-CPA SCHEME 
A. Main procedure 

As already mentioned, in order to find the statisti- 
cal averages lfTT|) and (|T2*|) we shall employ the single- 
site self-consistent approximation known as BEB CPA. 25 
This theory is capable to address simultaneously both 
diagonal and off-diagonal disorder, provided the second 
enters the random Hamiltonian in a generalized multi- 
plicative form. The last condition is intrinsically fulfilled 
in the excitonic problem under consideration because the 
coupling @ is given by a bilinear combination of the 
random p" with •& c ^ n being translationally-invariant de- 
terministic matrix. The BEB extension of the scalar 
CPAiS is based on the so-called BEB transformation. 25 
The latter represents a vector generalization of the mul- 
tiplicative ansatz originally implemented by Shiba^ to 
construct a single-site self-consistent theory for bond- 
disordered alloys with mean-geometric relationship be- 
tween hopping integrals. Upon such transformation the 



problem of finding the resolvent of a Hamiltonian with 
both diagonal and multiplicative off-diagonal disorder is 
equivalently reformulated for an operator which acts in a 
space with additional tensorial dimensions but contains 
only site-diagonal disorder. With application to Frenkel 
excitons the BEB transformation has been already real- 
ized in Sec. Ill Bl by considering instead of operators B\ 
(B n ) the polarizations P" (P,?^ ') equipped with vec- 
tor indices. As a consequence the dynamical equation for 
the scalar GF, 

G nm (z) = g n {z)S nm + g n {z) ^2 JniGi m (z), (23) 

i 

is replaced with that for the GT, 

^nm(z) = Jn(z)5 nm + 7n(z)'^2^nlTl m (z). (24) 

/ 

[Here and below omitting the indices we assume the usual 
rules for multiplication and inversion of tensorial quan- 
tities.] Unlike the case of Eq. (|23|l . the disorder enters 
Eq. I|24|) only in the site-diagonal form, namely, through 
the tensor 

lf{z)=p«g n {z)pi (25) 
associated with the bare local GF 

9n{z) = —^—. (26) 
z - e n 

At the same time, because the algebraic structure of 
equation l|24|) remains similar to that of (|23[1 . the usual 
CPA self-consistency arguments can be employed to ap- 
proximate the statistical averages defined in Eqs. 
and (|T2|) . The BEB-CPA scheme 2 ^ is accomplished via 
straightforward generalization of the scalar-CPA equa- 
tions within the tensorial formalism as outlined below. 

The disorder effect on the single-particle quantities is 
accounted in CPA by replacing the random local GTs 
"fn^{z) with a deterministic site-independent coherent- 
potential GT 7^(2). The disorde r-averaged GT in the 
momentum representation is found from Eq. (|24|) as 

f fe (z) = [l- 7 (z)^]- 1 7 W, (27) 

where 

< = ^E e " ifc(n " ro) ^ (28) 

is the momentum-space coupling kernel. In turn, the 
conditionally-averaged on-site GT T"P(z) is approxi- 
mated by assuming that the site, carrying its random pa- 
rameters, is placed into the same coherent environment as 
the one of Eq. (|27|l . Thus, finding this conditional aver- 
age reduces to solving a single-impurity problem realized 
by Eq. 124f) in which the bare GTs 7^ (z) for to ^ n are 
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replaced by 7 a/3 (z). A simple derivatiori22i2^2Si2i leads 
to 

f nn (z) = [1 - 7„(z)E(z)]~ 1 7 „(z), (29) 

where the site-diagonal self-energy T, a ^(z) describes 
coupling of the selected molecule with the coherent- 
background molecules constituting the rest of the chain. 
Clearly, being a functional of only the coherent-potential 
GT ^"^(z), this self-energy is the same as the one enter- 
ing the complete disorder-average of the local GT 



f nn (z) = [l- 7 (z)E(z)r 1 7 (z). 



(30) 



Provided both 7 a ^(z) and T"^(z) are represented by non- 
singular matrices one gets 



S(z)= 7 - 1 (z)-f^(^, 



(31) 



where the local GT is found in terms of 7 Q ^(z) from 
Eq. H27|) after the momentum-space integration: 



^nn (z) 



-[l-TtzRrSW. (32) 



The unknown coherent-potential GT ^f al3 (z) is to be 
determined in a self-consistent way. In spirit of the 
original CPA, one demands the single-site expectation 
of the conditionally-averaged GT to coincide with the 
completely-averaged local GT: 



(33) 



The latter guarantees that the conditional and the com- 
plete average of the GT found above will provide identical 
estimates for the single-site quantities such as LDOS or 
the local (single-molecule) polarizability. 

Summarizing, the BEB-CPA procedure amounts to 
solving equations which, including auxiliary definitions, 
are listed in formulas from (|25[) to J33J). For nontrivial 
disorder models and realistic forms of the coupling 
this can be done only with the use of numerical methods. 
An example of such numerical solution is given in Sec.lVl 
[see also Appendix . The analytic treatment of the 
low-disorder limit of the theory is presented in Sec. IIVI 



B. Analyticity and accuracy of BEB CPA 



As clear from the general discussion of Sec. Ill Bl the 
BEB CPA allows to find configurational averages of the 
exciton DOS and linear polarizability of a disordered 
molecular aggregate. It is natural to check the valid- 
ity of such nonperturbative theory asking, in particular, 
how close are the estimates to the actual quantities and 
whether they meet certain fundamental physical require- 
ments. This question can be partly answered by analyz- 
ing the structure of the BEB-CPA equations as done in 
a number of publicationsi22*2£i2§i2L2£ For completeness, 



let us outline the important facts concerning analyticity 
and accuracy of the BEB CPA. 

Gonis and Garland^! have proved that T"^(z) found 
within BEB CPA possesses the same properties as a func- 
tion of the complex energy z as would have the disorder- 
averaged GT calculated exactly. Namely, it is analytic in 
the whole plane excluding branch cuts on the semi-axis 
Im z = 0, Re z > 0, while the tensor 



1 

2?ri 



ff"(*) 



(34) 



that at z = lo + i0 + provides the exciton spectral density, 
is positively (negatively) defined at Im z > (Im z < 
0) . As a consequence the disorder-averaged polarizability 
will preserve causality, while the DOS and absorption will 
be nonnegative. 

The authors of Ref. |2f| have studied the accuracy of 
the theory to reproduce the energy-domain moments of 
the spectral density, each given by a coefficient in the 
high-energy expansion of the GT: 



-a/3 



^ z s+r 

s=0 



r af3 



dujLj s Af(Lj). (35) 



It was shown that the BEB CPA provides correctly the 
first three moments of the spectral density, 



Lkfl — lo, £fc,i — h + lo$klo, 
Lk,2 = h + h"dkh + k"dkh + lo'&klo'&klo, 



(36) 



where the corresponding moments entering the high- 
energy expansion of the local coherent-potential GT 
7°^ (z) are given by 



7,d 



77 r\1c 

(37) 

The fact that the number of these moments is high 
enough guarantees fulfillment of the sum rules for the 
DOS and polarizability, correctly accounting for redis- 
tribution of the total oscillator strength between com- 
ponents of the polarizability tensor in the presence of 
orientational disorder. It follows also that using BEB 
CPA one gets an exact value of the absorption lincwidth, 
provided the latter is defined from the weighted sec- 
ond centered moment of the spectral density [rj a ^ ~ 

\j L li ,2 L tt 0,0 - ( L fe£o,l) 2 /£fe£ ,o]- This meanS tliat ' 

even though there exists no unique way to define the 
width of an inhomogcneously-broadened absorption line, 
the BEB CPA still gives a reasonable estimate for this 
characteristics of the spectrum. 
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n-m>0 



IV. 



LOW-DISORDER LIMIT 



n-m<0 



FIG. 1: Schematic picture of the integration contour in the 
complex momentum plane. Each couple of vertical lines with 
opposite orientations give zero contribution to the integral. 
The crosses mark poles of the momentum-space GT. 



For a disorder-free chain the GT obtained from the 
BEB CPA coincides with the one known exactly. In order 
to establish the relation of our approach with the existing 
scaling theories of disordered one-dimensional excitons 5 
let us solve the BEB-CPA equations in the limit of weak 
disorder. A similar calculation for the bond-disordered 
binary alloy has been already done in Ref. |28j. However, 
due to the rather involved structure of BEB CPA the gen- 
eral treatment of the low-disorder regime is still missing. 
For the sake of simplicity we shall also relax generality 
imposing certain symmetries under which the theory is 
formulated in terms of diagonal tensors. 



Let us also mention that the theory does not seem to 
violate the fundamental inequalities \E— e n \ < J2 m \Jnm\ 
to be fulfilled for every eigenvalue E. These, in particu- 
lar, impose the lower, min(e„ — J2 m l«^nm|)j ancl the up- 
per, max(e„ -I- J2 m l^nm|)> bounds of the spectral region, 
where minimum and maximum are taken over all disorder 
realizations. Even though no rigorous proof of this prop- 
erty is known for the case of BEB CPA, it was shown to 
hold in the numerical solution of the BEB-CPA equations 
for some random alloy models. 28 As will be demonstrated 
in Sec. El these spectral bounds are not violated in the 
case of orientational disorder either. 

Another important aspect of the theory is the behavior 
of the disorder-averaged GT in the complex momentum 
domain. Considering (z) for complex k provides addi- 
tional understanding about the exciton coherence length. 
Since Tu (z) remains the same as Re A: is shifted by in- 
teger of 2ir, the integration contour in Eq. (|20|l can be 
modified as shown in Fig. ^ Provided the function is 
meromorphic (i.e., has only isolated poles) in k, the in- 
tegration reduces to sum over poles in the upper (lower) 
half- plane for n — m > (n — m < Q). [Due to the mir- 
ror symmetry of the problem each pole has a counterpart 
located symmetrically with respect to the origin k = 0.] 
To establish connection with the coherence length defined 
above it is sufficient to notice that only a single pair of 
poles will contribute to the asymptotics 121|) . namely, 
having minimal distance from the real axis Imfc = 0. 
Applying these arguments to the disorder-averaged GT 
in the BEB-CPA form (|27|l one concludes that the com- 
plex wave numbers £ Q/3 (z) of Eq. lj2"Tl) are to be selected 
among the solutions k of the characteristic equation 



det [1 - 7 (z)i? fe ] = 0. 



(38) 



A. General 

To proceed one needs first to specify the "low" disorder 
in terms of its probability distribution, entering the BEB 
CPA through the self-consistency equation (|33|) . We do 
this imposing smallness of the centered moments of the 
probability density which characterize deviations of the 
on-site energies and transition dipoles from their average 
values e — (e n ), p a — (p"}- The conditionally-averaged 
local GT in Eq. (|3*3l can be represented as a power series 
in e„ — e and — p a , so that its nth-site expectation can 
be written as 



(39) 



A vth term in the last sum is proportional to the vth- 
order centered moments of the probability density. The 
first three terms are given explicitly by 



(Tnn(z)) = f nn (z) 



(f nn (z)) 1 =0 
(f„n(*)) 2 = \ 



d 2 f nn (z) d 2 f nn (z) 



(40) 



a, (3 



de 2 

n ap d 2 r nn (z) 
' 2 dpfM 



where in the last equation 



a-2 = ((e„ - e) 2 ) , 



((e n -e)(i£-p«)> 



(41) 



Here we assumed that the Fourier-transformed kernel 
•d^ is analytic in k to provide meromorphic T^(z). This 
assumption is valid, for instance, if the coupling extends 
only over a finite number of sites. However, the last equa- 
tion may not lead to a correct coherence length for a more 
general nonanalytic coupling. 



are the second-order centered moments. The terms not 
included in Eq. I|4U|I will be neglected following our as- 
sumption about weakness of disorder. Later it will 
be demonstrated that the information contained in the 
third- and higher-order centered moments is indeed not 
preserved in the low-disorder limit of BEB CPA. 
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We shall additionally impose some simplifying restric- 
tions on the form of disorder and intermolecular coupling 
to be considered below. Namely, let the system to be 
such that (i) in some (orthogonal) basis of the three- 
dimensional vector space the coupling fl^m ^ s diagonal 
in the upper indices for all pairs of molecules, and (ii) 
the disorder probability density is symmetric under re- 
flections of two basis vectors. We shall label tensorial 
components corresponding to either these two basis vec- 
tors by symbols "J-i", "J-a") an d those of the remaining 
(third) vector by "||". The coupling energies J nm and, 
hence, the scalar GF G nrn (z) remain invariant, while ev- 
ery off-diagonal tensorial component of the GT T^(z) 
changes its sign under one of the above reflections [be- 
cause from a ^ (3 it follows that either a or (3 is different 
from ||]. Therefore, configurational averaging with a sym- 
metric probability distribution will result to a diagonal 

FS&(*)- 

Regarding the BEB-CPA theory, from Eqs. (J37JI and 
(|3U[I it follows consecutively that both the coherent- 
potential GT -f a/3 (z) and the self-energy T, af3 (z) reduce 
to the diagonal form. As concerns the low-disorder limit 
to be studied here, with our symmetry requirements the 
average dipole moment p a is directed strictly along the 
|j-axis, the vectorial parameter b% defined in Eq. I|41|) has 
nonzero component only along the same axis, while the 
tensor is diagonal. Using the explicit dependence of 
the conditionally-averaged GT on the nth site random 
parameters [see Eqs. I|25[) . (|26[1 and (|29p] the nontrivial 
(diagonal) components of the tensors (|40|l can be found 



fL(*))„ = 



o z - e - £ll(z)|p| 5 



r£i(*)> =0, (42) 



[z-e-Z\\(z)\p\*Y 
[z-e-£ll(z)bP]- 



c\(z-e)[z-e + iY,HzM 
[z-e-£H(z)H 2 ] 3 
4^ U (z)\p\ 2 



(43) 



y 



2 " 



Combined with the BEB-CPA self-consistency equation 
[Eq. (|33|) ] these formulas provide the starting point for 
the analytical treatment of the low-disorder regime. 

As seen from Eq. 1)43(1 . the solution behaves quite dif- 
ferently in the two distinct spectral regions. Namely, 
if z is sufficiently far away from the singularities gen- 
erated by denominator z — e — S"(z)|p| 2 one can look 
for the unknown j al3 (z) in an approximate form, replac- 
ing in JHI the BEB CPA self-energy S^fz) with its 



expression £ Q/3 (°)(z) for the disorder-free system. The 
solution found in this way will be valid only for energies 
not very close to the band edges of the disorder-free sys- 
tem. Around these points one can no more neglect the 
disorder-induced corrections to the self-energy compo- 
nent ElK°)(z) in Eq. lj4"3|) . Nevertheless, the last situation 
is still addressable analytically because near the band 
edges an explicit asymptotic relation between T"P(z) and 
7 a/3 (z) can be used. 



B. Iterative solution 

Away from the band edges of the disorder-free chain 
we can look for the coherent-potential GT in the iterative 
form 



*fP{z) = 7 (0)Q/3 (z)+7 (1)Q/3 (z). 



The first term with components 
7 (o)ll^ ) = J^i 



7 (0)±< (*)=0, 



(44) 



(45) 



corresponds to the disorder-free system and j^ a/3 (z) is 
a small correction to be determined. From equation 1)30(1 
linearized in the parameters 1(41(1 one has 



7 (i%) = 



£<°>ll(z)|p| s 



7 



{z-ef 



(46) 



where Tnn*^ {z) denotes the local GT on a single impurity 
site in the disorder-free chain, related to the disorder-free 
self-energy ^^(z) by Eq. ffSjl. 

The computation of the disorder-free self-energy 
YS°> a P(z) is not straightforward due to the fact that com- 
ponents 7 < -°^ i (z) and f„„ ' (z) vanish leading to an un- 
determinate expression in Eq. (|30|l . To overcome this dif- 
ficulty one can employ the following arguments. Consider 
the disorder-free chain with a single impurity molecule 
placed on the nth site in such a way that its dipole mo- 
ment p" has no component along the ||-axis. Since in 
this case the transfer energy J nm = J2 a Pn^nmP^ van ~ 
ishes, the nth site, decoupled from the rest chain, will 
be characterized by the local GT Tnn ^(z) equal to the 
bare quantity p"g n (z)p^. From Eq. J2HJ), for any of such 
p°, it follows that E a K s(0)Q/3 (^)Pn = 0. Therefore, be- 
ing a diagonal tensor, the self-energy will have nonzero 
projection only on the ||-axis, orthogonal to the subspace 
of the considered p". Using Eq. I|31|l for the nontrivial 
component of the self-energy we end up with a simple 
expression 



1 

w 



Gnn (z) 



£( )^(z)=0, (47) 
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where 



dk 



2ir 



(48) 



denotes the on-site scalar GF of the disorder-free chain. 
The first of Eqs. (|47|l represents a natural relationship 
between the scalar self-energy z — e - l/G { °l(z) and the 
main component of the tensorial self-energy in the ab- 
sence of disorder. 

Having derived the disorder-free self-energy we get the 
desired correction H46|) to the coherent-potential GT in 
the form 



7 d)x i(z) 



a 2 \p\ 2 G { V(z) 26>| 



4Gl°2(z) 



z). 



{z-ef 
3 



that usually represents an extremum of the bare exciton 
band. Hence, the iterative solution found in this section 
does not describe the most important part of the absorp- 
tion profile for polarization directed along the average 
dipolc moment. 



C. Scaling solution 

Let us now analyze the BEB-CPA equations near the 
lowest band edge of the disorder-free system. We shall 
use the standard (for J, but not H aggregates) assump- 
tion that this point corresponds to the center k = of 
the Brilloiun zone, thus coinciding with the energy ( de- 
fined in (|50|l . Around z — C it is convenient to represent 
the 1 1 -component of the coherent-potential GT as 



(49) 



\p\ 2 



a(zY 



(51) 



The final expression for the GT f ^ (z), valid in the low- 
disorder regime away from the band edges, is obtained 
by inserting the total j a(3 (z) of Eq. (JSJ) into (|57f) . 

Let us comment on the analytical properties of the 
obtained solution. The correction l|49|l obeys the large- 
z asymptotics 7 



Assuming the absolute value of the new unknown func- 
tion a(z) to be small compared to the total width of the 
exciton band one can keep only a few significant terms 
in the expansion of the coupling kernel d\ around k — 
while performing the momentum integration in Eq. 13211 . 
Provided the constant 



(i) a p 



(z) ~ c^/z. Hence, the disorder- 
(p a p /3 + c% )/z. Using 
definition (|41ll for it is easy to see that the coeffi- 



averaged GT behaves as 1^ (2) 

a/3 



j = d 2 4 



dk 2 



l/f 



(52) 



k=0 



cient p a pP + £f coincides with that of the first term has a finite (positive) value, that is the case if the real- 



(L^q = {PnPn)) m the general high-energy expansion 
B35J1 . That is, even though for the zero-order solution 
the asymptotics f 1 ^(2) ~ p a p 13 / z is not accurate, the 
correct form is restored as the first-order iterative cor- 
rection is taken into account. Furthermore, the next 
two coefficients (£2^, L^f 2 ) are also reproduced in the 
found solution provided one neglects the third- and the 
higher-order centered moments of the disorder distribu- 
tion function in Eqs. (|37|) . Nevertheless, the found GT 
can have in general spurious poles outside the real axis 
around the band-edge singularities of the disorder-free 
GF [Eq. (|48|l ] . Such nonphysical behavior is related to 
the fact that around the band edge the second term in 
the expansion l|44|l is no more small as compared to the 
first. For the two components f k i (z) the portion of spec- 
tral weight taken up by the spurious poles vanishes upon 
decreasing the disorder strength c 2 4 ■ This guarantees the 
correct reproduction of the optical absorption spectra for 
the corresponding polarizations, as well as the DOS away 
from the band edges. The same is true for the component 
r|(z) for a generic nonzero momentum k. On the other 

hand, the spectral density component ^5j!(z) at momen- 
tum k = is itself strongly concentrated near the energy 



space coupling $„ m decays faster than \n- 



at large 



intermolecular separation, we can use the effective mass 
approximation for the disorder-free exciton dispersion 



(53) 



This leads to the asymptotic expression for the corre- 
sponding component of the local GT valid around the 
band edge: 



y/2Jy/t-z + a(z) 



(54) 



Furthermore, in accordance with Eq. I|31(l . the respective 
self-energy can be represented as 

1 



(50) 



£ll(z) = — z- e - o(z) - V2J^C -z + a(z) . (55) 

\pV l j 

Concerning the two components T^(z), even though 

they share the singularity present in the above f j! n (2i), 
one can demonstrate that it is still possible to use the 
simple relationship f^(z) = 7 ±; (z), or equivalently to 
sctE ± '(z)=0. 

Writing the CPA self-consistency equation one can 
retain only the most singular terms of Eq. I|43|l near 
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the band edge z w £ where for vanishing disorder also 
o{z) 0. The equation for 7 reads simply 



7 ^(z) = 



(56) 



As for the component 7" (2), the only essential are 
the first and the third terms in the right-hand side of 
Eq. . By straightforward algebra we arrive at the 
self-consistency condition 



a(z), (57) 



2jy< - z + <7(z) 



that represents an equation to be solved for a(z). In this 
formula we have used the shorthand notation 



«2+44(Co)>l 



(58) 



for a combination which can be thought of as an effective 
disorder strength parameter governing the low-disorder 
limit of BEB CPA near the band edge. 

The algebraic structure of equation lfST|) is still too 
complicated to obtain its solution in a closed form. At the 
same time the scaling of o{z) with the disorder strength 
A can be found easily. Namely, in the interval \z — (j < 
J(A/J) 4 / 3 the solution (both real and imaginary part) 
behaves as 



<j{z) ~ J 



(59) 



Employing this scaling relation one can find the low- 
disorder asymptotic expressions for the absorption 
linewidth and the exciton coherence length in the res- 
onance region. Substituting (|51|l into Eq. H27|) the 
main component of the disorder-averaged GT around the 
phase-space point z w £, k rs can be written as 



z-C-(t{z) - Jk 2 /2' 
The corresponding polarizability 

\P\ 2 



xHz) = -M 



z-C-a{z)' 



(60) 



(61) 



obtained from the basic relation l|19|) , is characterized by 
a resonant behavior near the band edge. The width of 
the resonance, estimated as 77" ~ — Imer(£), scales with 
the disorder strength A as 



A 



4/3 



(62) 



Furthermore, one can make use of the effective mass ap- 
proximation (|53|l to get a scaling of the exciton coher- 
ence length in the spectral region around the absorp- 
tion resonance. According to the arguments of Sec. lIIIBl 



applied to expression l|60|l . the complex wave number 
k = £"(z) governing the asymptotic behavior of the 
disorder-averaged GT at large intermolecular separation 
is found from equation 



z-(- a{z) 



J, 



0. 



(63) 



Taking zw(wc end up with the following scaling of the 
coherence length in the resonance region: 



« ~ '7 



-2/3 



(64) 



It should be also mentioned that combining Eq. (|62|1 with 
(|64|l leads to the scaling relation 



77" 



J(7V" 



(65) 



which does not involve the disorder strength A. This 
dependence is universal in the sense that it follows only 
from the shape of the exciton dispersion in the disorder- 
free system in the vicinity of the band edge. Indeed, the 
scaling {55j l can be established without rigorous compu- 
tation of the spectral density but by simply noticing that 
77" and 1/N" can be considered, respectively, as uncer- 
tainties of the exciton energy and momentum brought 
about by the disorder scattering. The relation between 
these quantities imposed by the dispersion i|53|) is of the 
quadratic form (|65|l . 

It can be noticed that the terms of expansion 139f) pro- 
portional to the higher-order centered moments of the 
disorder distribution will bring stronger singularities in 
the vicinities of the band edges compared to ones of the 
second-order centered moments. Nevertheless, the dom- 
inant contribution to the solution will come solely from 
the second-order moments. Let us illustrate this for the 
simple situation when only diagonal disorder is present. 
A derivation analogous to that leading to Eq. I|59|) shows 
that the centered moment a„ of order v > 2 alone would 
produce the scaling <r v (z) ~ J(a„/ J") 2 /( 1 +^). If the 
probability to encounter a random energy e n falls off 
rapidly enough away from the mean value e we have a v ~ 
a 2 /2 . Thus, the contribution a v {z) - J(a 2 / J 2 )^/( 1 +") 
will be always suppressed by that of the second-order 
moment a{z) ~ J(a 2 /J 2 ) 2 / 3 considered so far. 

When applied to a molecular chain with only diagonal 
disorder the asymptotic expressions l|62l) and l|64(l coin- 
cide with those obtained in Ref. using different scaling 
arguments. Furthermore, for the case of diagonal disor- 
der the asymptotics (|62|) has been found in Ref. |22j by 
means of the scalar CPA. Our theory thus provides an 
extension of these results treating both the on-site en- 
ergy disorder and the disorder in transition dipoles on 
an equal footing through the combined parameter A of 
Eq. . The fact that BEB CPA is able to reproduce the 
known scaling relations is not surprising because this self- 
consistent approximation becomes exact in the disorder- 
free limit. 
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<|) (in-plane rotation) 



(chain's axis) g 

e ± =e 2 

FIG. 2: The studied geometry of orientationally disordered 
molecular chain. The dipoles are distributed in the plane built 
on vectors ejf, e± and forming the angle 8 with the chain's 
direction ef. 



V. ILLUSTRATIVE APPLICATIONS 

A. Model and simplified parametrization 

In this section we demonstrate how the general the- 
ory discussed so far can be applied to a specific model of 
disordered aggregate: a chain with purely orientational 
disorder in the transition dipoles. In addition to having 
constant absolute values the random dipoles will be as- 
sumed to lie in parallel planes forming an angle 9 with 
the chain's direction [Fig. [2] . The random in- plane ori- 
entation angles <f) n will be modelled by a box probability 
with density 1/(2$) in the interval \<j> n \ < $ and zero 
outside. We shall consider the molecules to interact as 
dipoles in an isotropic background, so that the tensorial 
kernel fl^m acquires the form 



i) 



li 



i) 



22 



v n 



i) 



:s:i 



-2V n 



12 



13 



2:\ 



0. 



(66) 



For the long-range (LR) dipole-dipole coupling one would 
have V nm = \n — m\~ 3 . However, since in this case the 
momentum integral (|32|) can not be reduced to a closed 
analytic form, we shall restrict ourselves to the more 
easily treatable nearest-neighbor (NN) [V nm = #| ra -m|,i] 
or next-nearest-neighbor (NNN) [V nm = £i n _ m M + 
^ln-ml,2/8] couplings. For concreteness we shall always 
assume that cos 2 9 > 1/3, so that the disorder-free chain 
($ = 0) represents a J aggregate. The singular case 
cos 2 9 — 1/3 is considered separately in Appendix [H] 
From now on the energy will be measured in the units of 
the NN interaction strength, and the dipoles will be cho- 
sen to have unit absolute values. We also shift the refer- 
ence point of energy to the position of the bare molecular 
level. 

Since the random dipoles are restricted in a single 
plane the theory can be conveniently reformulated in 
terms of tensors projected onto the corresponding two- 
dimensional subspace. We choose the new basis as 



= sin f 



ei = 0, 



= 0, 



= 1, 



cost 



e'i = 0, 



(67) 



i.e., the first vector is parallel to dipoles of the disorder- 
free chain, the second is orthogonal to it and to the 
chain's axis. In terms of the projected components 
Pn = cos<t> n , = sinr/>„ the three-dimensional dipoles 
are given by 



Pn 



ejfpjl + elpi 



(68) 



The advantages of this representation come from the fact 
that our disorder distribution density is invariant under 
the change of sign of p„. Hence, since upon projection 
the coupling l|66|) remains diagonal, 

tfL = (1 - 3 cos 2 B)V nm , ^ m = V nm , (69) 

the reduced BEB CPA will involve only diagonal ten- 
sors. As a result we end up with the affordable equations, 
whose explicit form for the case of NN and NNN coupling 
is given in Appendix^ Moreover, as the additional sym- 
metry requirements formulated in Sec. IIVI are fulfilled, 
the low-disorder solution derived there can be employed. 
Since all components of a physically-measurable tensor 
can be straightforwardly reconstructed using Eq. I|68|) 
from the two diagonal components of the correspond- 
ing projected tensor, the last two will be presented every 
time to illustrate the outcome of the theory. 



B. Numerical results 

Due to their explicit form the equations of Appendix lAl 
can be efficiently solved using standard numerical meth- 
ods. In order to test the accuracy of the theory we 
also performed exact diagonalization for an open chain 
of J\f = 250 molecules with either NN or LR coupling. 
The statistical error was reasonably small after averaging 
over 5000 disorder realizations with the energy domain 
divided into 400 output intervals. 

The tensorial components of the coherent-potential GT 
are plotted against energy in Fig. This figure can 
be used to illustrate analyticity of the theory and to 
predict certain features of the physical quantities to be 
presented below. Specifically, the functions lml/^(uj) 
and Iml/7 (o;) obey the correct sign that guarantees 
the positivity of DOS and optical absorption. From 
Eq. (|27fl it follows that the resonant denominators of 
the exciton spectral density components are given by 
[Re 1/7(0; ) — $fc] 2 + [Im l/7(w)] 2 , where polarization la- 
bels || or _L are assumed. Hence, provided Im 1/7(0;) <C 1, 
the spectral density will be characterized by a sharp res- 
onance located around the frequency u> to be found from 
Re 1/7(0;) — dk = 0. The solution of this equation for ei- 
ther polarization and momentum k can be found graphi- 
cally from Fig.[3Jc,d) [where, in particular, the positions 
of d\ and tf£ at k = are shown with horizontal lines]. 
Actually, the assumption on smallness of the imaginary 
part of the inverse coherent-potential GT is valid only 
in the case of component 7^(0; ) and only for moderately 
weak disorder. As concerns 7 _L (o;), the above condition is 
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FIG. 3: The imaginary (a,b), and real (c,d) parts of the BEB- 
CPA self-energy for the case of NN coupling. Several values 
of the disorder parameter are used: $ = 0.6 (solid lines), 
$ = 0.9 (dashed lines), $ = 1.3 (dotted lines), $ = 3.14 
(dash-dotted lines). 8 = 0.26 for all curves. The horizontal 
lines in panels (c) and (d) mark the diagonal elements of the 
NNN coupling matrix at momentum k — 0: 1?|_ = —3.603, 



never satisfied, so that one ends up with no well-defined 
quasiparticle peak in the corresponding component of the 
spectral density. 

The profiles of the disorder-averaged exciton DOS are 
presented in Fig. 21 Except for the well-known band- 
center Dyson singularity 3 ^ in the case of NN coupling and 
a similar but nonsingular feature for the LR coupling 34 
the numerical DOS is reproduced well within BEB CPA. 
As a confirmation of the analyticity of the theory no un- 
physical behavior of the DOS has been observed. The 
theory also captures properly the exact symmetry of the 
spectrum present in any tight-binding system without di- 
agonal disorder and with only NN coupling. The asym- 
metric DOS in the case of LR interactions is satisfac- 
torily approximated already by that of the NNN cou- 
pling. Furthermore, it can be observed that the upper 
and lower estimates of the spectrum region mentioned in 
sec. nun are not violated in the presented solution. For 
instance, in the case of NN coupling the spectrum should 
not extend beyond the interval \u>\ < 2max | J nn+ \\. The 
maximal absolute value of coupling accessible by varying 
parameters cj> n within the interval \<f> n \ < $ is given dif- 
ferently in the distinct domains of the parametric space: 
if cos 2 9 > 2/3 it equals to 3 cos 2 6—1, while in the op- 
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FIG. 4: The disorder- averaged exciton DOS for several values 
of disorder parameter $ calculated with BEB CPA (dashed 
lines) or exact diagonalization (solid lines), (a) Case of NN 
coupling (both BEB CPA and exact diagonalization). (b) 
Case of NNN coupling (BEB CPA) or LR coupling (exact 
diagonalization). 9 = 0.26 for all curves. The arrows mark 
the spectral bound estimates for the case of NN coupling and 
$ = 0.9, namely, u = ±3.60. 



posite case it becomes |1 — 3 cos 2 9\ cos 2 $ -I- sin 2 $ for 
$ < tt/2, or 1 for $ > tt/2. 

The disorder-averaged absorption spectra (normalized 
by the number of molecules in the chain) along the two 
essential polarizations, et? and e" , are presented in Fig.0 

The component Imx" [oj), corresponding to the direction 
of preferable orientation of the dipoles, demonstrates be- 
havior typical for disordered J aggregates: a sharp res- 
onance acquiring inhomogeneous width upon increasing 
disorder is located near the bottom of the excitonic band. 
As for the polarization orthogonal to the preferred orien- 
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FIG. 5: The excitonic absorption spectra for NN coupling 
and several values of disorder parameter $ calculated with 
BEB CPA (dashed lines) or exact diagonalization (solid lines). 
9 — 0.26 for all curves. The arrows mark estimates of the 
spectral center for $ = 0.9: ft 11 = -2.77, Q 1 - = 0.46. For 
<E> = 0.6 the absorption spectra resulting from the low-disorder 
limit of BEB CPA are presented (dotted lines). 



tation of the dipoles, the absorption component Im^ (w) 
is generated solely by orientational disorder. In contrast 
to the previous case no narrow resonance is observed, but 
rather a broad spectrum shifted to energies higher with 
respect to the bare molecular level. The positions of the 
absorption maxima can be roughly estimated from the 
zero and first moments of the spectral density as fi" ~ 

4=0,1/4=0,0' ^ ~ 4=o,i/4=o,o- Usin S E ^ © for 
these parameters in the simplest case of NN coupling 

we end up with QH = (1 - 3cos 2 6>)[l + sin2$/(2$)], 
ft 1 - = 1 — sin 2<5>/(2<£>). These expressions one more time 
confirm that the centra of the absorption spectra are lo- 
cated at essentially different energies for the two polar- 
izations. 

To provide additional understanding of the features 
observed in the absorption spectra, the complete BEB- 
CPA spectral density is plotted in Fig. EJa,b). For the 
considered disorder parameter (<3> = 0.9) the component 

aJJ.(oj) is characterized by a well-pronounced quasipar- 
ticle structure with most of the spectral wight located 
around the disorder-free exciton branch to = As for 
the part A^(u>), similarly to the already considered case 
of momentum k = (the only directly accessible in the 
linear optics), the function has no resonant behavior in 
all the phase space. The spectral density transformed 
to the real-space representation is shown in Fig. HJc,d). 

The component A| m (w) is essentially nonzero in the in- 
terval |n — m\ < 100, and for fixed exciton energy has 
oscillatory behavior as a function of intermolecular dis- 
tance n — m. Starting from zero near the bottom of the 
band the number of nodes monotonically increases as one 
moves towards the upper band edge. Since absorption is 
obtained from the real-space spectral density after the 
summation (XEJ, the absence of oscillations explains the 
strong absorbance for the | |-polarization near the band 
bottom. As concerns polarization orthogonal to the pre- 
ferred orientation of the transition dipoles, the function 
A^ m (u>) is confined on a few (|n — m\ < l-j-2) sites. This 
implies the absence of the coherent optical response of the 
molecules resulting to a weak and very broad absorption 
for this polarization direction. 

Figure illustrates variation of the exciton coher- 
ence length across the spectrum. The BEB-CPA co- 
herence length has been found extracting the parame- 
ters £, a/3 (z) of the large-intersite-separation asymptotics 
(|21|l from Eq. (|38[) . Owing to the symmetry of the prob- 
lem the latter splits into the two independent equations 
l — 7ll(2)i?| = 0, 1 — 7^(2;)^^ = 0. Consequently, each of 

the two components f \ m {z) and r^ m (z) is characterized 
by its own coherence length, to be denoted respectively as 
N^(ui) and N^iui). From the results presented in Fig.0 
it follows that N^(u>) strongly depends on energy, having 
maximum around the band center and decreasing as uj 
approaches the edges, whereas N^~{u)) is almost constant 
across the spectrum. Beyond the spectral region the pa- 
rameters £ a P(uj) are imaginary, and hence, the coherence 
length is finite, even in the uniform chain. Remarkably, 
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FIG. 6: (a,b) The tensorial components of the BEB-CPA exciton spectral density for the case of NN coupling. Only positive 
momenta k are shown. (c,d) The disorder- averaged spectral density in the real-space representation for a few fixed energies u. 
The BEB-CPA results (dashed lines) are confronted with ones of the exact diagonalization (solid lines) . The parameters are 
<5> = 0.9, 9 = 0.26. [The strong discrepancy at uj = is the effect of the Dyson singularity not reproduced by the BEB CPA.] 



upon variation of the parameter $ the two lengths N" (cj) 
and N (to) evolve differently: the first decreases while 
the second grows with disorder. Regarding the exact di- 
agonalization approach, the coherence length can be ex- 
tracted directly from envelope of the disorder-averaged 
real-space spectral density in accordance with Eq. (|2T1 . 
Unlike the case of DOS and absorption, much larger sta- 
tistical error is present owing to the fact that the quantity 
to be averaged (i.e., a tensorial component of the spec- 
tral density) depends on both energy and the real-space 
variables. Here one should keep small the width of the 
energy integration interval [5ui sa 1/ 400 of the total band- 
width in our simulations] to ensure that the additional 
decoherence due to interference of waves separated by 



energies smaller than Suj does not exceed the disorder- 
induced effect. If v(uj) is the excitation group velocity 
the last condition can be expressed as v(lo)/8u; » N^(lj). 
This inequality is not easy to fulfill around the band cen- 
ter [iVH(a;) is large] and around the band edges [v(uj) is 
small]. Apart from mentioned the numerical procedure 
suffers from the finite-size quantization taking place as 
the coherence length approaches the total length of the 
chain. 

For the model under consideration we have also ver- 
ified the accuracy of the low-disorder analytic solutions 
found in Sec. IIVI Since no diagonal disorder is present 
the centered moments and 6f of the distribution van- 
ish, and the only parameters which control the strength 
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use the iterative solution derived in Se c. HVBl w hich re- 
sults to expression Imx J -(o;) w Mc^ / v'C 2 — u)2 - I n con ~ 
trast, the component Imx"(k>), with most of the oscil- 
lator strength concentrated around the lower band edge, 
should be calculated using the approach of Sec. IIV CI In 
particular, taking into account that the disorder param- 
eter is given by A ~ $ 2 , the scalings (JH2Jl and 
readjyll ~ $ 8 / 3 , AT 1 1 ~ $- 4 / 3 ; respectively. The linewidth 
77" and the exciton coherence length at the absorption 
maximum extracted from the numerical solution of the 
BEB-CPA equations are plotted against $ in Fig. 02a, b). 
The presented profiles confirm the fulfillment of the de- 
rived scaling laws. Moreover, as illustrates Fig. [51c), the 
universal scaling (|65|l is observed in the whole range of 77" 
and TV" accessible by varying $ from to 7r, even though 
the two above dependencies hold only for <!> sufficiently 
close to zero. 



C. Experimental relevance 
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FIG. 7: Energy dependence of the exciton coherence length 
for several values of the disorder parameter <£>. The BEB-CPA 
results are presented in the case of both iV" and N ± (dashed 
lines) and ones of the exact diagonalization (solid lines) only 
for AT". 6 = 0.26 for all curves. 



of disorder are c\ ~ $ /45 and qr ~ $ 2 /3. To evalu- 
ate these constants we have used Eq. 11411 and retained 
only the leading terms in the limit of small $. In the 
magnitude of the averaged dipole moment one has to 
keep also the next-to-the-leading term, that results to 
the estimate \p\ w 1 — $ 2 /6. In Fig. [5] (panels with 
$ = 0.6) the analytic low-disorder absorption spectra are 
confronted with those calculated via direct solution of the 
BEB-CPA equations and by exact diagonalization. For 
polarization perpendicular to the preferred orientation of 
dipoles the spectral weight of the band-edge singularities 
(located at ±< = ±2(1-3 cos 2 9) for NN coupling) is neg- 
ligible compared to the integrated spectral weight. To 
approximate such wide-bandwidth absorbance one can 



In a typical experiment one deals not with a single 
molecular chain but rather with a macroscopic ensemble 
of such objects. The treatment reported here is appli- 
cable only if the coupling between molecules of distinct 
aggregates is weaker than the disorder-induced broad- 
ening in a single chain. The broadening mechanisms 
not related to the considered types of disorder can be 
accounted for via a phenomenological on-site dephasing 
rate 77, e.g., due to exciton- phonon scattering. The re- 
sults of the BEB-CPA theory can be straightforwardly 
extended to include this case by adding imaginary contri- 
bution it] to the complex energy variable z in all relevant 
functions. 

In most of the samples accessible in experiment, such 
as frozen solutions containing molecular aggregates, be- 
cause of the randomness in the chain's orientations one 
is not able to measure the individual spatial components 
of the microscopic (single- chain) polarizability tensor. 
Nevertheless, the available types of surface-deposited 
films 35 with strictly in-plane alignment or recently syn- 
thesized bulk materials with smectic-like 3 ^ alignment of 
aggregates still allow for independent observation of cer- 
tain spatial components of the single-chain polarizability. 
With a polarization-resolved spectroscopic study of such 
systems^ one can probe the orientational disorder within 
a single aggregate. In particular, a material obtained by 
stacking the chains described in Fig. |3 preserving per- 
fectly parallel alignment of the axes and planes represents 
a biaxial medium. The optical axes are given by the di- 
rections ejj 1 , e", with the corresponding polarizabilities 

proportional to x"^) an( i X^i^)? an d the one orthogo- 
nal to them, along which the optical response vanishes. 
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FIG. 8: Dependencies of the absorption linewidth (a) and the 
exciton coherence length at the resonance (b) on the disorder 
parameter <E>. (c) The same data is plotted in the coordinates 
"absorption linewidth — coherence length" . The parameter 9 
is 0.26 (solid line), 0.52 (dashed line), 0.65 (dotted line), 0.78 
(dash-dotted line). 



VI. CONCLUSIONS 

In summary, using BEB CPA we have investigated the 
linear optical properties of a disordered molecular ag- 
gregate with random single-molecule excitation energies 
and transition dipoles. The advantages of the used self- 
consistent approximation are its analyticity, exactness in 
the low-disorder limit and the capability to provide not 
only the DOS and absorption spectra but also the coher- 
ence length. To the best of our knowledge, the last quan- 
tity has never been addressed in excitonic systems with 
the coherent-potential approximation. Whereas the DOS 
is on-site characteristics while the absorption spectrum is 



related to the spectral density at momentum k = 0, the 
calculation of coherence length requires knowledge of the 
spectral density at a generic k. That is, the last quantity 
carries essentially new information on the single-particle 
properties of the system with respect to that contained 
in the first two. 

The interesting feature of the considered problem is the 
tensorial structure of equations brought about by orien- 
tational disorder in the transition dipoles. Such type of 
disorder leads to redistribution of absorbance between 
different polarization components with quite distinct be- 
havior of the corresponding coherence lengths. 
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APPENDIX A: BEB-CPA EQUATIONS FOR 
ORIENTATIONAL DISORDER WITH NN AND 
NNN COUPLING 

In the case of box probability distribution the config- 
urational averaging in the self-consistency equation ilMMI) 
can be performed analytically leading to 



(EH - S x )f L = q\ 



^eT 



(2 J --S")ri, = W iL ^r-'- 



where function q(z) is to be found from 
tan(g$) 



'z-E J 



tan<f>. 



(Ala) 
(Alb) 

(A2) 



From Eq. Q31[) the projected components of the self- 
energy can be expressed as 



E" = i 



f II 

J- nn 



i - r^- (A3) 

/ nn 



Performing momentum integration in Eq. (|32|l one re- 
lates the nontrivial components of the disorder-averaged 
on-site GT to the coherent-potential GT as 

FL = 7 ll /((l-3co S 2 0) 7 U) ) f4 = 7 ± /( 7 x ). (A4) 
The function f(x) entering these formulas is given by 

1 



Vl - 4x 2 ' 
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and 



V4 + 9x 



(V8x + V4 + 9z) 2 - 2x 
(V&x - \fl + 9xj 2 - 2x 



-1/2 



for the cases of NN and NNN coupling, respectively. 

For every disorder parameter $ and tilt 9 the equation 
presented in this Appendix can be solved numerically for 
7"(z) and r y ± (z). The example of such solution is shown 
in Fig. El 



APPENDIX B: FLAT-BAND SINGULAR 
SPECTRA 

In this appendix we study the single-particle proper- 
ties of the system described in Sec. IV Al when the an- 
gle 6 formed by the plane of the transition dipoles with 
the chain's direction is such that cos 2 9 = 1/3 (i.e., 
6 = 0.955...). The disorder-free situation for this value 
of 9 can be considered as an intermediate case between 
that of J and H aggregate, because the transfer energy 
© generated by the dipolar coupling l|oT)|) vanishes. Ne- 
glecting higher multipole contributions to the transfer 
energy as well as all other broadening mechanisms which 
go beyond the present model, let us concentrate on the ef- 
fect of orientational disorder to produce a finite excitonic 
bandwidth. [Considering instead the diagonal disorder 
alone results to a trivial physics: the off-diagonal part 
of the Hamiltonian is always zero and the molecules re- 
main uncoupled.] Our aim is to find using BEB CPA 
the asymptotic form of the GT around z = 0, where 
the DOS and certain components of polarizability tensor 
have a singularity. [Such singularity due to specific ar- 
rangement of coupling and disorder should not to be con- 
fused with the Dyson singularitj^, not accessible within 
BEB CPA] . For simplicity we shall consider only the case 
of NN coupling. 

From Eq. HA4jl with cos 2 9 = 1/3 it follows that 

?lL(z) = 7" (z). According to Eqs. i|A3(l the corre- 
sponding component of the BEB-CPA self-energy is zero, 
= 0. On the other hand ^ ± (z) is nonzero at any 
finite z. Taking this into account it is easy to conclude 
that the auxiliary function defined by formula ljA2f) be- 
haves as q(z) — > 7r/(2$) upon z — > 0. Applying the same 
limit to Eqs. (|Alb|l we get 



1 



1 

2' 



7 (z) ~ - — W— . 



(Bl) 



(B2) 



Combining formula i|Bl|) with Eq. (|Ala|l we obtain the 
asymptotics of the remaining unknown functions at z — > 
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FIG. 9: (a) Exciton DOS and (b), (c) absorption spectra in 
the case cos 2 9 = 1/3. The BEB-CPA results (dashed lines) 
are confronted with ones of the exact diagonalization (solid 
lines) and the derived asymptotic expressions (dotted lines). 
The disorder parameter "1> is equal to 0.9 for all curves. 



I, in 2i 
Az) = tHz) <~ \ — . 



(B3) 



It turns out that, irrespectively on the strength $ of 
orientational disorder, the component f|l„(z) of the 
disorder-averaged on-site GT is singular as z — > 0, i.e., 
at the energy of the bare molecular transition. Strictly 
speaking, upon introducing some degree of diagonal dis- 
order or a finite inelastic dephasing rate rj such singular- 
ity would be suppressed. Nevertheless, one can still apply 
the derived asymptotic formulas if the absolute value of 
z is much larger than 77 but smaller than the total band- 
width generated by orientational disorder. 

The exciton DOS, related to the on-site GT found 
above by formula i|17fl . shows an analogous singularity 
around uj = 0: 



1 1 



4<I> 



(B4) 



As concerns the absorption spectra, for the two essential 
polarization directions these are characterized by com- 
pletely distinct behavior: 

Im X l! H~^^i Imx^M ~ (B5) 



The DOS and the absorption profiles calculated via nu- 
merical solution of the BEB-CPA equations and ones 
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of exact diagonalization are confronted with the derived 
asymptotic formulas in Fig. EI AH approaches agree well 
close to the band center. Noticeably, unlike the usual 
case cos 2 9 ^ 1/3, the absorption component Imx(w) is 
symmetric with respect to the central energy u> = (still 
only the NN coupling is assumed). 

Regarding the exciton coherence length, since #| = 
it follows that AT II (z) = for every z. In other words, for 
polarization along the preferred orientation of the dipoles 
the disorder-averaged spectral density depends on the 
intermolecular separation as Imf| m (z) ~ S nm , i.e., we 
encounter completely incoherent optical response of the 
molecules. For the orthogonal polarization in the vicinity 
of z — we have 




(B6) 



N ± (lu) 



4$ 1 



(B7) 



As a result the coherence length N (w) is finite within 
the band but diverges at its center. In the vicinity of lu = 
the corresponding real-space spectral density behaves 
as Imr„ TO (0) ~ cos(7r|?i — to|/2), i.e., coincides with the 
center-of-the-band spectral density of a finite-bandwidth 
uniform chain. 
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